QUANTUM DIFFUSION, MEASUREMENT AND FILTERING 

V P BELAVKIN 

Abstract. A brief presentation of the basic concepts in quantum probability 
theory is given in comparison to the classical one. The notion of quantum 
white noise, its explicit representation in Fock space, and necessary results of 
noncommutative stochastic analysis and integration are outlined. 

Algebraic differential equations that unify the quantum non Markovian dif- 
fusion with continuous non demolition observation are derived. A stochastic 
equation of quantum diffusion filtering generalising the classical Markov filter- 
ing equation to the quantum flows over arbitrary *-algcbra is obtained. 

A Gaussian quantum diffusion with one dimensional continuous observation 
is considered. The a posteriori quantum state difusion in this case is reduced 
to a linear quantum stochastic filter equation of Kalman-Bucy type and to the 
operator Riccati equation for quantum correlations. An example of continu- 
ous nondemolition observation of the coordinate of a free quantum particle is 
considered, describing a continuous collase to the stationary solution of the 
linear quantum filtering problem found in the paper. 
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1. Introduction 

Beginning in the mid seventies, modern probability theory has, (along with tra- 
ditional subjects, such as dynamical systems with random perturbations), been also 
concerned with fundamentally new stochastic objects — quantum dynamic systems 
with an inherently probabilistic nature. Mathematically, the concept of quantum 
probability arises not because of the lack of information for a complete description 
of the object, the instability of chaotic motion or the inaccuracy of measurement 
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but is due to the noncommutativity of the algebra of random variables which are 
represented by the operators in the Hilbert space. As quantum probability theory 
is an intrinsically stochastic theory, it cannot be stated within the framework of the 
Kolmogorov axioms [J , which assume the fundamentally deterministic description 
of the classical systems under the given point states to G fl. It is based on differ- 
ent axioms [2], [Hj namely, the Neumann axioms, whose greater generality can be 
demostrated even in the case of a finite number of alternative elementary events 
u> = 1, . . . , n. 

Let us illustrate for this simple case how the classical probability space (fi, J 7 , P) 
can be represented as a special case of the quantum space that is defined by the 
triple (H,„4, E). Here H is a (finite-dimensional) complex space of column- vectors 

h = [rf], rf e C with scalar product (g\h) = £?W = 9*h for g = [C] G H 
defined by the weights (probabilities) pi > 0, A is an associative, but not necessarily 
commutative matrix algebra X — [£ l k ] closed under the involution X i— > X* defined 
by the Hermitian conjugation 

(X*g\h) = {g\Xh),[e k ]* = \pT 1 ^iPk], 

and with matrix / = [S\,] as the identity I £ A, and ELY] — (e\Xe) is a positive 
normilized functional (ELY*A] > 0, E[J] = 1) of expectation of noncommuting 
variables X defined by a fixed unit vector e 6 H, where ||e|| 2 := (e|e) = 1. 

Classical random variables x : f2 — * C can also be described by the multi- 
plication operators X = x , (xh)(u>) = x(aj)h(oj) in the complex Hilbert space 
H = L 2 (£l, J 7 , P) of T -measurable P-square-integrable functions h : £1 — > C , 

\\h\\ 2 = J \h{uj)\ 2 P{duj) = (h\h) < oo. 

Their expectations E[S] = J x(to) P(du>) are defined as (e\xe) by the unit function 
e(ui) = 1 which is normalized with respect to any probability measure P. Thus, 
however, only commutative operator algebras A are obtained, whose elements are 
given in the finite case of — {1, . . . , n} by all the diagonal matrices x — [£,{i)S l k ] 
with the commutative product, corresponding to the pointwise multiplication of 
the functions x(u>) — £(i), where lu = i. 

Conversely, any quantum probability space (H, A, E) can be reduced to the clas- 
sical one (fi, J 7 , P) only in the case of the commutativity of the algebra A; in the 
finite-dimensional case this is realised by simultaneous reduction of the commutting 
matrices X € A to diagonal form [£(*)<5fe] • The probabilities pi of the elementary 
events u> — i in the diagonal representation are defined by the restriction pi = E[Pj] 
of the functional ELY] = ^ £ {i)Pi on projective matrices Pj = (i) 5*] , ^ (i) = $) . 

In this article a quantum analog of diffusion and the problem of its continuous 
measurement and stochastic filtering, that gives the solution of the Zeno paradox 0] 
(as a result of establishing an a posteriori stationary state), are considered within 
the framework of the noncommutative algebraic approach. A derivation of the 
stochastic equation is given for a nonnormalized a posteriori quantum state, which 
is obtained in [5] by renormalising the basic equation of nonlinear quantum filtering 
[Hj. The solution of the equation has been found for the case of linear quantum 
diffusion of canonical commutation relations, obtained previously for the quantum 
Gaussian case by means of linear Markov filtering methods in [7] , |S] • 
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In presenting the second (basic) section, we deliberately avoided the questions 
concerning the sufficient conditions for the dense definition of the unbounded infini- 
tesimal generators that guarantee the uniqueness of solutions of quantum stochastic 
and operator equations; this is beyond the scope of this article. We only point out 
that in the first and second sections these questions are not relevant (see [5]) for 
the Markovian case with complete pre-Hilbert domain D in the initial Hilbert space 
H, corresponding to boundedness of the operators L and H in D = H. Moreover, 
a solution exists for an unbounded algebra A of canonical commutation relations, 
which is considered in the third section, in the framework of a quantum calculus of 
kernels for the operators L, H £ A in the initial Fock scale {Fj|^ > ll^U] if their 
inductive limit UF^ is chosen as D. Besides, the explicitly solvable model of this 
section with linear unbounded generators L and H , does not require the estimates 
obtained in these scales. 

For completeness the notation and explicite methods of quantum stochastic in- 
tegration and the proof of their estimates in Fock scale [5] are briefly presented in 
the Appendix . The comprehensive statement of the author's general approach, 
outlined above, and the estimates for the integrals can be found in JHj, [2 • The 
earlier results on quantum stochastic calculus in the framework of Hudson and 
Parthasarathy approach |12| . are reviewed in|13|. 

The approach presented generalizes the results for purely quantum diffusion in 
[T^| to the case of an arbitrary initial algebra A. This enables a unified description of 
quantum and classical diffusion, their observation and filtering as special algebraic 
cases. In the sections 3 and 4, a one-dimensional variant of an infinite-dimensional 
quantum Gaussian filtering j!4j is presented as well as an example of observation 
of a coordinate of a free quantum Brownian particle; this was analysed earlier in 
|16j by the method of solving the a posteriori Shrodinger equation |15|. 



1.1. Basic Notation. Let H be a complex Hilbert space and D C H be a dense 
subspace defined as an inductive limit (see appendix 1) of some scale {H^S, > 1} 
in the space H. Let the initial algebra A of noncommutative random variables 
describing a 'quantum object' at the initial moment t = be represented by an 
involutive subalgebra A C B(D) of linear operators X : D — * D, X* £ A, having 
(an inductively) continuous conjugate X* : D — > D, (X*x\^j) — (x\Xip) with 
respect to the scalar product in H, with an identity operator I £ A. 

Let us denote by H — H<g> F the tensor product H and the Fock space T — T(K.) 
over the Hilbert space K, — L 2 (R+) of a 'quantum noise' Wt(g), g £ K, and let the 
pre-Hilbert space T> be an inductive limit of the Hilbert scale = H^ <g> , £ > 1 , 
where is the natural Fock scale (see appendix 2) over fC. We shall consider 

the quantum noise as a set {w t {g)\g £ L 2 (R + )} of Brownian motions t i— > Wt{g), 
represented in T by self-adjoint operators 



with a Gaussian state on the algebra generated by them, which is induced by the 
vacuum function 8$ £ T . Here {a ri a*\r £ R+} are canonical operators of creation 
a* and annihilation a r in J- (see Appendix 3) called quantum stochastic integrators, 
and g*(t) = g~(t) . Note that each operator function 1 1— > Wt{g) that has commutative 
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values [w s (g), wt(g)] — is equivalent to a classical Brownian motion with intensity 
|g(i)| 2 , with respect to the vacuum vector e = 5$. This follows from the formula 

e iw(g) _ e ia*(g) e -i||g|| 2 e ia(g*) 

and e a ^^0 = 5$ for any / G K. , due to which the quantum characteristic function 

E [ e «M9)] = («5 |e™*W<J ) 
coincides with the clasical Gaussian characteristic function 

|2, 



/ 6XP ^/ 5 , ( r ') dw r}P(dw) = exp J 



\g(r)\Mr 



of the standard Wiener process w t . However, the different Brownian motions w t (f) 
and w t (g) with f*g ^ g* f do not have any classical representation on a single 
probability space (f2, J 7 , P) because of noncommutativity (see A3 in Appendix 3): 

[w t (f), Mg)} = J q (j{r)g(r) - g(r)f(r))dr. 

Definition 1. Let {At\t £ R+} be an increasing set of involutive subalgebras 
At C A s , t < s of the operators X t e B(V) generated by operators X £ A, w t (g), 
g £ K,, such that 

X t eAt*> [X t) Y] = 0, Vy e B(I>) : [X , Y} = 0= [W t {g), Y] , 

where the operators Xq £ Ao, Wt(g), g S L 2 (0, t] are assumed to act in H = H x T 
as Ao ® 1 and I ® w t {g). A measurable operator function F(t) : T> — * T> is called 
adapted if -F(i) S ^4* for almost all t £ R+. 

We shall consider here only quantum stochastic integrals of the form 

(2.1) i\{F,D)= I (F(r) dA r + D{r) dA*), 



where A r — I®a r , A* — I® a* and F, D are locally square-integrable (see Appendix 
4) together with adjoint F*,D* adapted operator-functions R+ — > At- Note that 
on the exponential vectors, described by the product-functions h(r) = /c®(t)?/>, 
where tp £ D and k £ L 2 (R + ), the itegrals (|2.1|l are weakly defined as the usual 
operator integrals 

r t / 

dr. 



0| il{F,D)h) = y [F{r)k(r)+D(r)k(r)]h 

This gives in particular, E[« (F, D)] = for E[A] = (e|Ae), where e(r) = ^(t)^ 
is the exponential vector, corresponding to fe = 0. Such operator integrals on the 
exponential domain were constructed by Hudson and Parthasarathy |12 | for the 
case of bounded F (t) ,D(t) . 

For the adapted integrals (|2.1|l . the quantum Ito formula 20| _ EI can be ob- 
tained as in the case of l|A.5(l (see Appendix 5), corresponding to Xt(t) = for 
almost all t. This formula defines (see Appendix 5) the pointwise multiplication 

X(t)Y(t) = X(Q)Y(0)+ f d(XY)(r) 



of the adapted operator-functions 

X(t) = X(0)+ f dX(r), Y(t)=Y(0)+ [ dF(r) 



o Jo 
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in terms of the product dX(t) dY(t) = D(t)*F{t) dt of their stochastic differentials 
dX = D*dA + D dA* , dY = F*dA + F dA* : 

d{XY) = dXY + XdY + dXdY = D*Fdt+ (D*Y + XF*) dA 

(2.2) +(DY + XF)dA*. 

The classical Ito formula for the stochastic integrals 

Ii(f,u) = [ f{r,uj)dw r 
Jo 

with respect to the standard Wiener process w = {w t \t £ R.+} can be obtained from 
()2.2|) by the Segal one-to-one transformation w : Wt i— ► w±, where uit = at + a^ . The 
latter represents the adapted operator integrals i Q (f, f) — J Q f(r) dw r for the non- 
anticipated functionals fit) = f(t,w) of commuting selfadjoint operators w = {wt} 
with 

ll/*0ll?:= f\\Kr)5,f<^ 
Jo 

in the form of the Ito integrals uj\l\Q , /)] = Iq(J, u>), so that 

H(f,f)s 9 \\ 2 = \\fs 9 \\f = J ||^(/^)|| 2 p(dc), 

where P is the standard Wiener probability measure. 

1.2. Quantum diffusion. Quantum stochastic evolution in the open system {At} 
is described by an adapted family {i(t)\t S R+} of ^-representations : X 
X (t) of the initial algebra A C B(D) into B(T>) , i.e. of linear maps A — > At, with 
the properties: 

i(t, X*X) = i(t, X)*L(t, X), i(t, I) = I := I <g> i. 

It is called diffusion motion if the operator- valued functions t i— » X(t) have the 
quantum stochastic differentials of the form 

(2.3) dX(t)+C(t)dt = D*(t)dA t + D(t)dA;. 

Here C(t) = j(t, X) is an adapted operator-valued function, locally integrable 
(p = 1) for every X E A defined by the linear maps j(t) : A — > At, t £ R+ . 
D*(t) = 5*(t,X), D(t) = S(t,X) are adapted operator-valued functions, locally 
square-integrable for each X £ A, defined by the linear maps S*(t), 6(t) : A,—> At. 

Define the output process Y = {Y(t)\t € R+}, which is subject to measurement 
and described by a commutative family of (essentially) selfadjoint operators Y(t) — 
Y(t)* on V with the initial condition Y(0) = and stochastic differentials 

(2.4) dY(t) = G(t) dt + F* (t) dA t + F(t) dA* . 

Here G(t) = G(t)* £ At is essentially self-adjoint and locally-integrable (p = 1). 
F(t) £ At is locally square-integrable together with its conjugate: F*(t) = F(t)* ; 
G(t) and F(t) are adapted operator- valued functions of t G R+. 

Unlike the classical case, not every involutive subalgebra B of A , but only 
a central one B c A fl A' , defines the conditional expectations E[X|B] for any 
state vector e £ H as the positive projections A — > B which are compatible with 
E[X] = (e\Xe) such that E[E[X|B]] = E[X] for all X £ A. Hence, not every 
stochastic process described by the equation l|2.4|) . can be considered as an output 
process for quantum diffusion, defined by equation (|2.3|l . but only that for which 
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the posterior expectations of X(t) with respect to the observation Y(s), s < t, 
exist. 

Definition 2. A process Y{t) is called causal, or nondemolition with respect to 
the process X{t) if 

(2.5) [X(t), Y(s)] := X(t) Y(s) - Y(s) X(t) = 

for all t > s, s £ R+. 

The nondemolition condition together with the self-nondemolition of Y , i.e. with 
the commutativity [F(i),y(s)] = 0, Vt, s, is necessary and sufficient JT] for the 
existence of the conditional expectations 7r(t,X) = ~E[X(t)\Bt] for the operators 
X(t) — i(t,X) , with respect to the ^-algebras 

B t = {Y £ B{V) | [X,Y] = 0, Vie B{V) : [X,Y(s)] = 0, Vs < t} 

generated by the family {y(s)|s < t} and for every initial vector-function eel?, 
||e|| = 1. If the process Y with Y (0) =0 is nondemolition with respect to the 
coefficients C, D* , D of the equation l|2.3J) . then it is nondemolition with respect to 
the solution X , corresponding to any initial X (0) S A . This and other sufficient 
conditions of the next proposition obviously follow from the integral representation 

X(t) = X- f (C{r)dr-D*(r)dA r -D(r)dA*). 
Jo 

Proposition 1. The integrals X (t) of (|2.3J) with X £ A are ^-representations 
i(t) \ X i ► X(t) iff the linear maps 

7(t):Xi — *C(<), — ►!>*(*), — ► £>(*) 
satisfy the following differential conditions 

(i) 1 (t,x*) = 1 (t,x)*, s(t,x*) = s*(t,xy, vx£A, 

(ii) 7 (i, = X)* 7 (t, X) + 7 (t, JC)*i(t, X) - 6(t, X)*S(t, X), 

S(t, X*X) = i(t, X)*S(t, X) + 5(t, X*) i(t, X) = S*(t, X*X)\ 

(iii) 7(t,J) = 0, *(*,/) = = tf* (t, J), V<eR+. 

The process X(t) satisfies the condition l|2.5|l iff the stochastic derivations C(t),D* (t),D(t) 
also satisfy the condition (|2.5I) as X(t) with respect to the nondemolition process 
Y(t) for all X £ A, and the derivatives G, F*, F in (|2.4|l satisfy the differential 
nondemolition conditions 

[X(t),F*(t)]=0 = [F(t),X(t)], VteR +J 

(2.6) D*(t)F(t)-F*(t)D(t) - [G(t),X(t)]. 

PROOF. The stochastic differentials dX(t) = X(t + dt) — X(t) of the linear 
*-maps i(t) : X i— > X(t) are defined by the linear *-maps ~f(t), S*(t), 5(t) by virtue 
of linear independence of the fundamental differentials dt, dAt and dA^. The 
conditions (ii) are found by applying the Ito formula H2.2|) to X(t)*X(t): 

d(X(t)*X(t)) = dX(t)*dX(t)+dX(t)*X(t) + X{tydX(t) 

= [S(t, X)*5(t, X) - 7(i, X)*i(t, X) - i(t, X)*j(t, X)] d t 

+ di (S(X)*t(X) + l(X)*S*(X), S(X*)l(X)+l(X)*S{X)). 

By equating the stochastic derivatives of this differential and 

di{t,X*X) = dt t (S*(X*X),S(X*X)) ~-f(X*X)dt, 



QUANTUM DIFFUSION, MEASUREMENT AND FILTERING 7 

we obtain that 6* and 6 are the derivations of the algebra A, and —7 has the 
positive-definite dissipator 

l(X)*^(X) + j(X)*t(X) - ~f(X*X) = S(X)*6(X). 

The condition (iii) follows from d bit, I) = because of the independence of 7, 5*, (5. 
If Y(t) is a nondemolition process for X(t), then 

[<Ur(i),y( a )] = [X(t + dt),Y(s)] - [X(t),F(a)] =0 

with t > s ; hence the nondemolition for C, I?*, 13: 

[C(t), F(s)] = 0, [D*(t), Y(s)] = 0, [D(t), Y(s)] = 0, Vi > a, 

follows by commutativity of y(s) with the indepedent differentials di, d^4t and dAj? . 
Applying equation (|2.2I) to the differential of the commutator [X(t), Y(t)] = we 
obtain (taking into account the equality [dX(t), Y(t)} = 0): 

d [X(t), Y(t)] = [dX(t), dY(t)] + [dX(t),dY(t)] + [X(t), dY(t)] 

= (D*(t)F(t) - F*(t)D(t) + [X(t),G(t)]) dt 

+di t ([X,F*],[X,F}) =0, 

which yields the differential self- nondemolition conditions Ij2.6|l . Hence, all condi- 
tions of the proposition are necessary. 

1.3. The Markov case. The quantum diffusion l|2.3|) with coefficients 
C(t) = i(t,C t ), D*(t) = t(t,U t *), D(t) = i(t,D t ), 

corresponds to the Markov stochastic evolution (in strong sense) . Here Ct , D% , D t € 
A are defined by the structural maps 

lt :X^C u 5* t :X^D*, S t : X h- D u 

for which the conditions (i)-(iii) indicated above were obtained by Hudson and 
Evans in ^3]. The self-nondemolition conditions 1)2. 6|) give the restrictions for the 
coefficients G and F* , F in this case. We shall restrict ourselves to consideration of 
the standard case F(t) = Iq = F*(t) of the indirect measurement 

(2.7) y(t)= I t(t, G r )dr + I(g)wt 

Jo 

of the diffusion of a square-integrable initial process defined locally by G t £ A over 
the standard Wiener process Wt represented in Tt by the operators I®wt = A t + A^ . 
It is not hard to prove [5] that Y(r) = Vt(I®w r ) V t *,Vt > r E R + by the uniqueness 
of the stochastic operator equation 

dVt + \G{t)V t dt = iG{t)V t du t , V = I, 

o 

where 

i 1 
G(t) = b(t,G t ), u t = - (a^ - at) = ^w t (i). 

The above implies the local unitary equivalence of the processes At + A* t and Y(t), 
which is always the case for the locally norm-square-integrable operator-functions 
G t :H^H|. 
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Corollary 1. In the case under consideration the condition l|2.6(l completely defines 
the structure of the inner derivation X t— > [X, Gt] for St — 5* t : 

(2.8) 6 t (X)= ±[X,G t ]-a t (X), 

where at (X * ) = a t (X ) * , V X E A is a ^-derivation of the algebra A. The conditions 
(i)-(iii) here also define the structure of the maps j t : A — > A in the form 

7t(X) = l -\ t {X)-P t {X), 

(2.9) \ t {X) = ~[G u [G t ,X]]+G t a t (X) + a t (X)G t -a 2 (X), 

where a 2 (X) = a t (a t (X)), f3 t : A — > A is some *-derivation 

/3 t (X*X)=X*j3 t (X)+/3 t (X*)X, MXeA. 
Proof. By taking into account the differentiation property 
[G t ,X*X]=X*[G t ,X] + [G t ,X*]X, 
and similarly for at(X*X) we obtain 

Xt(X*X) = l[G u X*[G u X] + [G u X*]X] 

+G t (X*a t (X) + a t (X*) X) + (X*a t {X) + a t {X*) X) G t 
-(X*a 2 (X) + 2a t (X)*a t (X) + a 2 t (X)*X) 

= X*X t (X) + X t (X)*X -^[X,G t }*[X,G t ] 

+ [X, G t ]*a t {X) + a t (X)*[X, G t ] - 2a t (X)*a t {X) 
= X*X t (X) + X t (X)*X-2S t (X)*6 t (X). 
Hence 7° = ^Xt possesses the property (ii) of the map 7 t : 

7?(X*X) = X* 7 °(X) + ^{Xyx 6 t (X)*5 t (X), 

and the ^-property 7° (X*) = 7° (X)* , as does j t in (iii). From here the result that 
j3 t = 7° — 7 t is a ^-derivation follows. 

The maps 7 t are called generators for the Lindblad equation 

d^/dt + S o 7t = 0, 

where fioj(X) = fi(j(X)), which is an algebraic analog of the Kolmogorov equation. 
This is satisfied by the operators Xq = ^(X), X E A of the conditional expectation 
Hq(X) = (diji\X(t) <5$) , with respect to the vacuum function 5® € T defined by 

X t i/j= [t(t,X)h] (0), V/i = V<8><50, VeD- 

The differential conditions obtained for Markov diffusion are necessary for the 
existence of a unique solution X(t) of equation (|2.3fl for all 1(0) = I £ i and so, 
for the Lindblad equation. They are sufficient in the case ^1] of constant at , /3 t and 
Gt and the boundedness of the algebra A (for instance, with D = H). Moreover, 
as was proved in [TJ]|, the maps bit) : X 1— > X(t) are representations of A in At, 
satisfying the condition of nondemolition i|2.5[l for all X E A. This is also true 
under significantly more general conditions of local p-integrability over the norm of 
the operators Gt : H — > H (with p = 2) and the maps at (p = 2) and (3 t (p = 1) 
from A C B(R) to A. 
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In the case of inner derivations 

a t (X) = i[S t ,X], p t {X)=i\H t + j{S t G t + G t St),X 

(as is the case for the von- Neumann algebra A), a quantum Markov diffusion is 
defined by structure maps of the type 

5 t (X) = [X,L t ], S* t (X) = [L* t ,X], 
(2.10) 7t (X) = i(L*[L t ,X] + [X,Li]Lt)+i[X,H t ], 

where 

L t =~G t + iS t) L* t = ^G t -iS u H t =HZeA t . 

3. Quantum diffusion and filtering. 

2.1. The a posteriori dynamics. The quantum diffusion i|2.3[l under nondemo- 
lition measurement (|2.4|l . is described by the classical random variables Xt(w) = 
(X(t)) t {uj) of the conditional expectations (X(t)) t — E[X(t)|£? t ] , for the operators 
X(t) = i(t,X) , on the trajectories w e O of the process Y(t) , with respect to 
initial vector- valued function eo=?p ®Sy l ,ip ED. As was established for the first 
time in [3], 0, the random process Xt : ^ — » C , considered as a stochastic map 
X(t) i— > Xt(u>), satisfies the Ito filtering equation for quantum diffusion 

(3.1) d(t(t,X)) t + (~f(t,X)) t dt = (n(t,X) - (G(t)) t i(t,X)) t dY 

with respect to the stochastic map (• ) t : X(t) i— > Xt- In this equation Y(t) = 
Y{t) — J Q (G(rf) dr is an innovation martingale for the observed process 12.4fl . and 
ft(t) : A — > B(D) is a linear *-map, which in the case of F(t) = Iq = F*(t) is of the 
particularly simple form: 

(3.2) K (f, X) = i (G(t) X(t) + X(t) G(t)) - a(t, X), 

where a(t, X) is defined by the ^-derivation 8(t) + 5*(t) — —2a(t). Equation l|3.1|l . 
derived in [S] by means of the martingale methods of quantum nonlinear filtering, 
extends the basic equation (8.10), JSI of the optimal diffusion filtering to the case 
of the noncommutative operator algebras A. By complete analogy with the clas- 
sical case the quantum Ito formula was used with the innovation process and the 
representation theorem which requires the conditional expectations x t {uj) to 
exist. This requirement is met by the self-nondemolition condition (|2.5() for all s < t 
which is trivially sutisfied in the commutative case s > t. 

In the quantum Markov case of the indirect measurement (|2.7|) the conditional 
expectation Xt(u>) = n t (X,uj) of the operators X(t) can be found as in the classical 
case by solving an autonomous stochastic equation for the a posteriori state TTt(X) = 
(t(t,X))t- The latter is defined on the trajectories w as a linear stochastic positive 
normalized map oj(TTt{X)) — x t (uj) of the algebra A into C satisfying the condition 

xt(u) y{u) P* (dw) = ( eo | X(t) Ye ), VXeA. 

In this equation, Y € Bt is any bounded operator in the algebra of the observed Bt, 
y{ui) = oj{y) is the Segal transformation of the operator y in T that corresponds to 
the unitary-equivalent operator I®y = V t *YVt, and Pq is an induced (by the unitary 
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transformation) probability measure on the trajectories u>|[0,i) := {u>,-|f G [0,£)} , 
restricted to the interval [0, t) with respect to the initial vector-state 

yMPS(dw) = (V t *M,<8>*B)) I (I®y)V t *(%®6 )). 

2.2. The filtering equation. Let us sketch the essentials in the derivation 
of a stochastic Markov quantum filtering equation, obtained for the general out- 
put process in [2]. First, we shall prove that the vacuum conditional expectation 
lAg(X) — (5q>\TTg(t,X)5$) of the product n g (t,X) = t,(t,X)e g (t), where 

(3.3) e g (t) = exp j J g{r) &Y{r) - \ g{rfdr 

satisfies the linear evolution equation 

(3-4) ^ ti* g (X) + fig o lt (X) = n g (G t -X - a t (X)) g(t), 

where fi^(X) — X, and G- X = (GX + XG)/2. Let us asume the uniqueness 
of this solution, which is always true for locally p-integrable bounded (over the 
norm) maps G±, a t (p = 2) and [3 t (p — 1). We shall prove that H g (X) is the 
mathematical expectation of the product fl g (X) — /i (X) e* , of the stochastic 
operators //'(-X", a>) = Lo[fi {X)\ (which satisfy a quantum filtering equationt) with 
the exponentials 

ft r i i 

e*(w) = exp J \g(r)dw r ~ -g(rf^dr = u^), 

(which are defined with respect to the trajectories w : t > w t of the standard 
Wiener process Wt, t S R+). The above means that the output process Y(r), 
restricted by any t G R+ , is absolutely continuous with respect to the standard 
restricted process w — {w r \r < t}. This follows from the unitary equivalence 
Y(r) = V t Y r V t *,Vr < t and Y r = I ® w r , Wt G R+ , representing the out- 
put up to a time t with respect to the initial vector-function eo = i/'o ® ^0 an( i 
e = V t *eo correspondingly. The probability density p (w) = Pg(dw)/P(du;) , 
for the measurement of the trajectory {uv|r < t} , of the process Y on the in- 
terval [0, t) , is defined with respect to the standard Wiener probability measure 
P(dw) by the formula Pq(ui) — ip^(I,uj). Here <Pq(u)) = <p a o ^(ui) is the stochas- 
tic functional iPq(X,u>) — (t/>q|/x*(.X", lu) i/j ), which corresponds to the initial state 
<p (X) = (tp \Xip ) on the algebra A. Finally, we deduce a nonlinear equation for 
the a posteriori state n t (u) using the ordinary Ito formula and the normalization 
of the stochastic functional <Pq(oS), 

Theorem 1. Let the equation \3.4-\j has the unique solution X g — (i JX), cor- 
responding to the initial condition (ig(X) = X for each X £ A and g G /C 
Then it coincides with the vacuum expectation X g — (5^\TT g (t, X) 5$), where 
7T g (t,X) = i(t, X)e g {i) and is defined by the Wiener average 

(3.5) X\ = J exp | J g{r) dw r - X - girfd r J P(dw) 

over the continuous trajectories uj G fl. Here -X^(cj) = w[/2 ] is the Segal transfor- 
mation of the solution X 1 — fi{X) to the operator filtering equation 

(3.6) dtf(X) + o lt (X)dt = fi t (GfX- a t (X)) dvl t 
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with initial condition Jl Q (X) — X. In this case the linear stochastic equation 
\H.6\) also has a unique solution in the Ito sense, which defines almost everywhere 
(Pq(uj) ^ 0) for each ip Q G H the a posteriori state 

where the probability density /Og(w) = ( - o |/ t (aj) ip ) is given by the positive operator 
= n (J, uS), satisfying the martingale property 

I t {uj)V(doj\w r ) = r(w),Vr < t. 

Proof. First we find a quantum stochastic equation for X g (t) — ir g (t,X) using 
the Ito formula 

de g (t)=g(t)e g (t)dY(t), e ff (0) = 1, 

where 

dY(t) = G{t)dt + F(t)*dA t + dA* t F(t). 
We obtain according to (|2.2(l 

d(X(t)e g {t)) = dX(t)de g (t)+dX{t)e g (t)+X(t)de g (t) 
= [ g D*F-C + gXG){t)e g {t)dt 
+d il (D* + X F*gD + X Fg) e g (i) 

= n g (t, {G t X +IG t )^- a t (X) g(t) - 7t (X)) dt 

+n g (t, 5* t {X) + g(t) X) dA t + n g (t, 6 t (X) + g{t) X) dA* u 
where the explicit form (|2.8I) has been used for 

D* = S*(X), D = 8{X), C = 7(X), GX + D = G-X -a(X) 
with F = I = F*. 

Taking into account the martingale property of the quantum stochastic integral 
(|2.1(l with respect to the vacuum-vector 5$ S we find the equation 1)3. 4)1 for the 
operator 

f i t g (X):Tp^S g (X)Tp = {ir g (t,X)hm 

in H is implied by the action of ir g (t, X) = X(t) e g (t) on h = ip (g) $0 : 

d M *(X) = M *((G t -X- a t (X)g(t) - lt (X))dt, VI G A 

where j t : A — > -4 is defined in the form 1)2.10)1 . 

Now, if -^'(cj) = ^(X, lu) satisfies the stochastic equation 

dX\uj) + C 1 (uj) dt = D'ioj) dw t , X°(uj)=X, 

we can derive a differential for 

M*(X)=M t (X)e* I e*=exp|^ (g(r)dw r - i ff (r) 2 dr)|, 

by means of the Ito formula of classical stochastic calculus. 

dX t (w) + C t (uj)dt = D t (oj)dw u X°(u)=X. 
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Using de* = g(t)e g dw t we have 

d(l~*e*) = dX*de* + dX*e* +l*de* 

= (<?(*) D* - C 4 ) e*dt + {& + X* fl (t)) e*dw t , 

what can be written in the form of the stochastic equation 

dXl M + (G* (w) - 23* (w) 5 (t))dt = (23* (w) + X* («) 5 (t))d^ 

for (w) = w ^JC*e* J. Hence, the mathematical expectation (13.51) of Xg(uj) = 

X t (uj)e t g (uj) with respect to the Gaussian measure P of the standard Wiener process 
w satisfies the equation 

d M * (X) = (Dlg(t) - C\) dt, fJ? g (X)=X. 

Comparison of this equation with equation l|3.4[l gives the coefficients 

g* = J c'M4HP(d w ), D*=y D ( H<( w )P(d w ) 

in the form: 

C*= M *(7 t P0), D\ = ^ g (G t -X-a t (X)). 

Consequently G*(cj) = w(G'), 23*(u>) = cj(23*) are the coefficients 

£*=/%*(*)), 23* = ju t (G r X-a t (X)), 

that define an equation for X*(u>) = ^{X.lo) in the form of 13.6(1 . The solution 
2*(w) = /u*(2, u;) to this equation for the initial condition X = I defines a positive 
operator-valued diffusive process 2*(u;) = w(ju (2)) which satisfies to the martingale 
equation d2* (u;) = G 4 (w) du>t with the initial condition 2° (w) = 2 , where G* (w) = 
/i* (G,v) and the properties a t (2) = = 7 t (2) are substituted into i|3.6|) . Thus 
the Theorem 1 is proved. 

2.3. The classical case. The remark that follows provides an explanation why 
equation l|3.6[) is a noncommutative analog of the Zakai filtering equation. 
Remark 1. Let A be a commutative algebra equivalent to the space C°°(R d ) of 
infinitcly-diffcrentiable functions x : R d — ► C with the pointwise product, and the 
involution x*{z) — x{z). Then equation (|3.6|) is an operator representation of the 
Zakai equation 

(3-7) d[x t Zo + -A tt ,i o dt=(g t +V t )fi t zo dw t , V° Zo =6 Zo , 

for the nonnormalized a posteriori distribution n z (dz,Lo) of the Markov diffusion 
process z{i) described by the stochastic equation 

(3.8) dz + c t {z) d t — a t (z) dv t , z{0) = z , 

with the indirect measurement 

dy(t) =g t (z(t))dt + dw t 
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defined by the standard Wiener process wt — —vt- Here 

/ x(z)V t fJ.(dz) = - ^a k (z)x' k {z)fi(dz), 

^ ^ k=l 

« d 

x{z)A ti i{dz) = +2 / ^c t fc (z)4(z) M (dz) 
J fc=i 



E 

k,l=l 



^ a t {z) a t {z) x'u(z) (i(dz), 



with x' k = dkX, x'li = dkdix. The integral J //* (dz, oj) — defines the proba- 

bility density of the output process y(r, ui) on the interval < r < t with respect to 
the Wiener distribution P(dw) with given initial state z G R d , 8 Z0 (dz) = 1 with 
z G dz, <5 Z0 (d;z) = 0, and z ^ dz. 

Indeed, in the case of the commutative algebra A ~ G°°(R d ) , G t is the multi- 
plication operator by the given function gt(z) of the state z G R d of the Markov 
process z(t) . This process has the generator r y t (X)(z) — \T t x\(z), defined by the 
diffusion operator Tt on the measurable functions x : z i— > a;(z), a; G C°°(R !i ). The 
indirect measurement of gt(z) is given by the output process 

V{t) = / Sr(2(r))dr + «! t . 



Since [Gt,X] = 0, VX G A, 5t = —a t is a real derivation 

a t (X") (z) = a t (z)dx(z) := a t fe (z) d k x(z), d k = 
and 7 t = G*- at — (3 t — \a\, where Gr at(X)(z) = <? t (z)at(X)(z), the operation 

(3.9) T t = (a t g -b t )d- i (a t 5) 2 = c*3 fc - i aJafoQ 

is a standard generator of the diffusion process z(t), with = a l t (g t Sf — ^dia^) — b k . 

Note that the noise v t — —w t in the classical system appeared essentially the 
same as in the observation channel because it was represented in the Fock space of 
the Wiener process Wt- In order to represent a classical stochastic system in the 
same way with the noise Vt =/= which is independent of Wt, it is necessary to start 
from the Fock space T — T(/C) over K, = C m (8 L 2 (R+) with multiplicity m > 2, 
as is the case in |14| . 

2.4. The a posteriori equation. In the general case the filtering equation (|3.6|l 
defines the nonnormalized a posteriori state tp — tp o/i , which is the vector-state 

<Po(X) = (ip \Xip ) for all (p (X) = (^> |Jf^ ) in the case of inner derivations H2.10|) . 
The normalized a posteriori state TTt(X) = ip (X) /Zp Q (I) satisfies the (nonlinear) a 
posteriori equation 

(3.10) d9 t (X) +9 t o lt (X)dt = 9 t (K t (X) - n t (G t ) X)dw u 
with initial condition ttq(X) — ip (X) , where 

K t (X) = G t -X-a t (X),dwt=dwt-*t(G t )dt. 

The nonlinear stochastic equation H3.10|l is the Markov case of the general a pos- 
teriori diffusive equation i|3.1|) with innovating martingale dY represented by dw. 



14 



V P BELAVKIN 



It can be deduced from the linear one (|3.6|) for the nonnormalized state <p t (X) = 
ip (jj, Q (X)) by applying the classical Ito formula to the product tp t {X) = 'p 7r t (X) 
and noting that the positive martingale /Sq = <Pq{I) has the stochastic differential 
d/?Q = Zp (G t )dw t . Indeed, 

dfioiX) = ^d7T t (X)+dp*7f t (X)+p*d7ri(X) 

= tp (G t )Ttt °K t (X)dt+ (cp (G t )nt{X) + p 9t dw t 

-?o(^t °7tP0 + nt oK t (X)TT t (G t )) dt 

= <p i (GrX-a t (X))dw t -<? Q o' rt (X)dt, 

where Kt(X) = n t (X) — Tr t (Gt)X. Note that in the deduction of the equation the 
following relation was used: 

pl^t ° K t (X) %{Gt) = TTt ° K t (X) tp (Gt). 

4. Linear quantum diffusion with observation 
Let S be a symplectic {(-space, i.e. a complex space with the involution 

and skew-symmetric bilinear ft- form: s:HxS->C, 

s(r], rf) = -s(?7 B , 77), s(rjt,r))* = s(r?, rf). 

We denote by Re 3 the real space of the jj-invariant vectors 
rj = if € S, and assume that Re 2 is a Hilbert space with respect to the scalar 
product (£, rj) = (77, £), satisfying the inequality 

where £ 2 = (£,£), V 2 = (ViV)- 

A linear map R : E — * /3(D), satisfying the ((-property R(rj)* — R(rf) defines an 
operator representation of the canonical commutation relations if 

(4.1) [R( V ),R(^)] = R( v ) R{tf) - R{r^)R{r,) = \ s( V , 77") I. 

on a complex pre-Hilbert space D. It is called Gaussian with respect to a normalised 
vector V € D, ||^ || 2 = (ip \tp Q ) = 1 if 

(4.2) ty I e *«ty ) = e^«)-« 2 / 2 e, 9 (0, £ G Re S. 

Here $0(77) = (Vi^o) is the linear continuous (J— functional 
^oiv)* — floiv^) °f the mathematical expectation (^ |i£(?7) ^ ) = i?o( 7 ?), that is 
defined by some ?5o € Re H by means of the complexified bilinear form 

on S. The product (£,77) corresponds to the symmetric covariance 

Re (R(£) % I R( V ) iP ) - MO Mv) = (6 »?>> 77 e Re 5. 

The exponents e R ^' = X(£) = e lR ^ are defined as unitary operators, which form 
the Weyl family {X(£)\£ € Re S}, 

X(0 X(77) - e ls «<">X(£ + 77), Vfc 77 e Re H, 
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with the self-adjoint generators 

fl(OV> = -i^*(WU=o, CeReS. 

Such a representation can be realised in the Fock space H = F over the comple- 
tion K of the (quotient) space S with respect to the (semi) positive definite scalar 
product 

Indeed, the space F can be defined as the comletion of the (quotient) span D of 
the exponential vectors {77® | 77 € S} with respect to the scalar product 

00 

( £® I ®) = g I )» = exp(C I ). 

z — ' TV. 

n=0 

Let ^4* : S — * 0(D) be a linear map, defining the creation operators in D by the 
adjoints A* (77) — A(rf )* to the annihilation operators 

A(i) : r?® h-v (£ I 77)77®, V£eReS. 

From the canonical commutation relations 

[A(r)t),A*(r))]=(r]\r))I>0, VijeH, 

one can obtain the relations l|4. 1|) for the linear combinations 

R(rj) =Mv) I + 2UA(rj), 2^A(rj) = A* (77) + A{rf). 

This defines the Gaussian representation £ 1— >R(£) with respect to the vacuum 
vector tp = 0® in F, so that A(rf)ip = 0, V77. The Weyl operators X(£) are 
defined in F as 

(4.3) X(0 = O (0 e A * W) e AW) , V£ e Re 5. 

We obtain the representation (|4.2|) : #o(£) = (^0 1-^(0 V'o) ^ we take account 
the fact that e A( ^tp = ip . 

Let us denote by j : 77 1— * jr](— 77) a canonical bounded map from the Hilbert 
space 5 with respect to the norm 

\rj\ = (77*, 77) 1 / 2 = y/ (Re ?7) 2 + (Im t?) 2 , Im 77 s = Re 777, 

into the pre-Hilbert space S with respect to the (semi) norm ||?7|| = (77 | 77) 1 / 2 , 

\\v\\ 2 = (v\ V) = \v\ 2 + T^s (77,77") = |?7| 2 + s (Re77,Im?7) 

< |?7| 2 -I- |s(Rer7,Im77)| < |t7| 2 + 2 | Rerj\ \lmr)\ < 2\rj\ 2 . 

Then we can write (£|t?) — (£ gry), where g = 1 — |s, so that gij is the complex 
bounded functional t?(£) = (^\r)) = (£,g?7) on ReS which together with = 
(^iSv)* = (v\0 defines the Hermitian functional 

2 Re 1? = t9 + ^ = 2 Re?/ + slimy, 

where s : ReS — > ReS is a skew-symmetric operator (£,s?7) = s(i],£), \srj\ < 2\r]\. 

Let us consider a quantum diffusion for the operators R{t, 77) = n(t,R(r))) with 
continuous indirect sequential observation of the operators 
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Here Lt = i?(Ct); ^t — R(Ct) are defined by the weak locally square-integrable 
families {( t }, {Ct} 01 the elements £ t , Ct £ 3, i e R + so that 

/ £ r (i? J ,tf)dr < oo,V?9 € E,t G R+, 
Jo 

where 

(4.4) ^^)=^(Ct)^(d) = |(^,C t )| 2 - 

Taking into account the fact that any derivation of the Weyl algebra A is internal, 
we consider only the structural maps p.lOJI , given by the operators L t , L\ and by 
a Hamiltonian H t in the initial Fock space H = F . Moreover, we shall assume that 
the Hamiltonian is obtained by the normal ordering H t =: ht(R) : of a quadratic 
form ht (R) of the operators A + A* such that 

{% | H t ^ v )=v t ('&) + ^ t (^d), = tf + 2Re(gry). 

In this equation ^ = exp{— \(j)\rf) + A*(r])}ip are the normalised exponential 
vectors 

^ = e-M a /y, M 3 = (^|^)=l, V?7 € 3, 
generating H = F as the completion of the linear envelope {ip v }- {vt\t € R+} is a 
locally integrable family of the linear j}— forms vt ■ 3 — >■ C, and {wt|t € R+} is a 
locally integrable family of real symmetric forms 

w t (#,<!?) = w t (Retf,Rei?) +w t (Imi?,Inn?) = w t (i?,#). 

We shall consider the forms Vt,u>t to be continuous, so that 

where u* € Re 5 and u>t is a symmetric operator which is bounded with respect to 
the norm in 3. 

Proposition 2 Under the assumptions made above about the linearity of Lt, L* t 
in {R(r])} and quadraticity of H t the equation l|2.3|) for = R(t, rf), r\ 6 3 with 
D = [Jf, i], I?* = and C defined by the generator (|2.1UI) . is linear and 

autonomous with respect to the family \R(j])\r\ € 3}, 

(4.5) dR(t, rj) + R(t, in t s7j) dt = du f (77) + Iv t (sr]) dt, 
where K t is a complex bounded operator in 3 and 

dv t (rj) = s(i],Im( t )dw t + s (77, 2 Re C t ) du t = - dw t (s((, irj)), 

where ut = ^A^ . 

Indeed, for the operators H t with the quadratic Wick symbols ht(i}) = t>t(#) + 
|wj(i?,)J) the following commutation relations hold: 

i[H t ,R(r))] = v t (sr]) I + R(u t sri), VryeS. 
In addition, the equations (|4.1|l also give 

i[R(r,),L t ]= S (r l ,( t )I, i[L* t ,R(v)] =4(lv)I 

for L t = R(( t ), L* t = R(d). By substituting this into fFTfy we obtain C t = 
R(iKtsrj) — vt(sr)) I, where 

1 j. ■ 8 
K t = 2 "ft + It =£*-£*■ 
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Hereet is an operator in S, that defines the positive form (|3.4|l 

£ t (#,i?) = (#,e t tf>, (#,efi?), = £ t (tf,#), VtfeH. 
In the right-hand side of the equation (|2.3() we obtain the quantum noise 
D* t dA t + D t dA* = [R(Ct),R( V )] dA t + [R(r,),R(Q] dA* t 

(4.6) = i(s(r),e t )dA t -s(r 1 ,C t )dA* t )=dv t (v)- 

The following theorem establishes the existence and uniqueness of the solution 
of equation ((4.5(1 together with the integral of the locally square-integrable real 
function g(t) with respect to the stochastic differentials 

(4.7) dY(t) = R(t,( t + (i)dt + dW tl 

where Wt — At + A* t is the Fock representation of the Wiener process — an error 
of measurement L t + LI in Fock space T — T(/C). 

Theorem 2. Let the equations j4-5\ ), 14 f or the linear diffusion be defined by 
the functions v t ,io t and e t , which are locally integrable with respect to the norm in 
H, such that 

nt rt 

h| t (1) =/ |u r |dr<oo, IkI?^/ |K r |dr<oo, Vt, 



where \v t \ = \Jv\, \k±\ = sup{Kt($',i?)||i9 / |, < 1} , and 

/ft \ 1/2 

IC + C'|| 2) = f j (2RcC) 2 drj <oo, VteR + . 
T/ien i/iey /iGwe a unique solution defined in the Hilbert space F ® J 7 as 

(4.8) i?(i,r?)+ / S (r)dr(r)= / s(v r , Vr ) dr + R( v (t)) + t (/*, /). 

Jo Jo 

by the quantum stochastic integral i\ ) {f* 1 f) = a{ft) + a *{ft) with 

o(/ t *)= f f*{r,rj)dA r , a*(f t ) = [ f(r,r,)dA$. 



Here rj(t) = <Po (t,rj) — Vo> Vr — ^^{tiV)) r G [0, 4[ is the solution of the backward 
conjugate equation 

(4.9) -f] r + iK r sr) r = g(r) (( r + £»), 

wif/j i/ie boundary condition r\ t —r\. The complex functions ft, ft of r E R+ (me? 
t)£S eguaZ zero at r > t and are defined for r < t by the real function g(r) as 

f*{r,rj) = g(r)+is(^\t, V ),(i) = /(r,^)*. 

Proof First, we write the weak solution of the equation 1(4.9(1 in the standard 
Duhamel form 

r lr = <j>(t)r ] + [ g(s)cf> r (s)(t s + (l)ds, 



where 4> r {t)n = (p^\t,rf) is the solution of the equation 1)4.9(1 with zero right-hand 
side g = 0. The resolving operator <p r (t) exists as the chronologically ordered 
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exponential 

4> r (t) = ^ J " J ^t 1 s---K tn sdt 1 ---dt„. 



n -°r<ti<---<t n <t 



This comes from the estimate of the norm 

\ct> r (t)\ = sup||0 r (t) 77I : \n\ < l} 

< £l S l" /■■■ / |Kti|"-|«t»|d*l-"dtr. 



™ _0 0<ti<---<t„<t 

< exp I y |2/c r | dr 



which is finite as |s| < 2. Hence, the linear form (r] r , $) is uniquely defined on 5 9 1? 
for every r € [0,i) as a bounded functional with the estimate 



1/2 



\(r) r ,ti)\ < |<0 r (t),0>| + / | 5 (s)|<<Ms)2ReC s) tf>|ds 

< H + | 5 |f) K + c itf)Q* |0T (s)l? | 2 d ^ 

< (hi + Iffl^lC + C S |f } |0| exp {2 IkI^}. 

Now we can integrate the left-hand side of the equation (|4.8|l by parts, taking 
into account l|4.6f) and (|4.9I) : 

77) + / g(r) (i?(r,2ReC r )dr + 25RdA r ) 
Jo 

= R(t,r))+ (2$t{g(r)dAr} - R(r,r} r -iK r sr) r )dr) 
Jo 

= R(0,r] )+ I (2^t{g(r)dA r }+dR(r,r] r ) + R(r,iK r sr] r )dr) 
Jo 

= R(4f\Ui 1 ))+ ( (2^{g(r)+is(r 1r ,Ci)dA r } + v r ( S r, r )dr). 



Here di?(r, rj r ) is the quantum stochastic differential dR(r, rj)]^—^ , satisfying the 
equation l|4.5|l for t = r. This proves Theorem 2. 

Remark 2 The solution R(t, rj) of the equation l|4.5|l given by the integral (|4.8|) 
for 5 = preserves the commutation relations (|4.1|1 and satisfies the nondemolition 
principle 

[R(t, v ),Y g (t)]=o, v 77 e s, s a 2 (R+), 

with respect to the commutative (self-nondemolition) processes 

Yg(t) = I 9(r)dY(r), g € L 2 (R+). 
Jo 

Indeed, by using the quantum Ito formula we can obtain 

[di?(t,?7),di?(i,?7 J )] = [d« t (?7),dW)] = 7i (sW) dt, 

[dfl(t,»?),dy fl (t)] = [d« t (?7),dw t . 9 (t)] = lS (77,C + C B )ffWdt. 
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Hence, if [R(t, 77), R(t, 57*)] = [l/i)s(r),rft)I, then 

77), d#(t, 77 s )] = -K t (s7j,s?7»)d* 



and 

dfiJ^tj),^*,^)] 

d[J2(t,»7),y ff (t)] 



K t - K t + 7t) ( s? 7, W) di = 0. 

[dfl(t,»7),y B (t)] + [R(t,r}),dY g (t)] + [dR(t,r,),Y g (t)} 
= { [R(t, v ),R(t, C t + Cf)] + fofa, C t + C?) l} 9(t) dt = 0, 
if [R(t,rj), Y g (t)] = and, consequently, [dft(t,rj), F g (t)] = 0. 

Example. Let us consider the simplest nontrivial case of the space S = C 2 of 



column vectors 77 
conjugation 77" = 



'/ ;J 



, r) p ,ri q £ C with the involution 77 
) and nondegenerate symplectic form 



77" (a complex 



s(77, 77 tt ) = 2(77^77* - 77 ? 77p. 

This corresponds to the canonical commutation relations [P,Q] = for the 

operators momentum P and coordinate Q for the one-dimensional quantum parti- 
cle. They are defined in the Fock space over K = C as (14.111 for R(rj) = r\ p P + rj q Q 
due to the degeneracy of the semi positive definite scalar product 



« I V) 



VpCp + VqCq + 2 S ( 7 ?' £ ) = (C + *fg) fap ~ "7?)- 

This Fock representation in F= T(C) is associated with the Gaussian state (|4.2|) 
(which corresponds to the standard scalar product (£",77) = £*£ p + ) and is 
equivalent to the Shrodinger representation in H = L 2 (R). In this representation 
Qo = x and P = f^, and 

^° {X) = (2^/4 6XP { ~ 2 {X ~ ^ + 2~ 

where q = (V> I Wo) = $ q and p = (t/) \Pip ) = "d q . 

Suppose that this particle exhibits free quantum Brownian motion, i.e. H = 
(1/(2to))P 2 is its Hamiltonian, where to > is the mass of the particle, so that 



p(t) = p®i+d t , 



1 



— / P(r)dr, 
m Jo 



Q(t) =Q® 1 

% = iy/\(A* t -A t )=I® 2VA3A t * 

is a realisation in the Fock space F eg) JF, where .F = r(L 2 (R + )), of the Wiener pro- 
cess with intensity 2A with respect to the vacuum-vector <5$ G T . This corresponds 
to the quantum stochastic equation 14. 5|) with the parameters Vt — 0, 

1/2to \ _ / -2 \ _ \/A / 

J ' St ~ ^ 2 J ' Ct "2 v l 

In this case the equation 1)4. 6fl describes an indirect observation of the coordinate 
of the particle, realised by measuring the increment of the commutative process 

Y(t)= ( y/\Q(r)dr + W u Y(0) = 0, 
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with precision A > 0, where Wt — At + A* t is the standard stochastic error of the 
measurement represented by another realisation in Y®T of the standard Wiener 
process with respect to 5$. Note, that by the noncommutativity 

[v(s), w(r)] = 2iV\ min{s, r} I <S> 1 

it is impossible to represent these quantum processes together as Wiener classical 
processes on the same probability space, although each of them separately allows 
such a representation. In addition, the nondemolition condition of the pair (P, Q) 
with respect to the observation Y can be verified directly by calculating the com- 
mutators 

[P(t),Y{s)] =0, [Q(t),Y{s)] =0, Vs<t. 
(One can check, that for s > t the above do not commute) 

5. Markovian filtering of Gaussian quantum process 

The solution (|4.8II , obtained for the linear quantum diffusion equation 1)4. 5fl with 
continuous observation l|4.6|l . enables us to find the Weyl operators l|4.3l) in the 
Heisenberg picture X(t,£) = L(t,X(£)). To this end let us represent the product 

X g (t, = exp | J (g(r) dY(r) - \ g{rfdr) J X(t, 

as the exponent of the operator 1)4.8)1 : 

X g (t,0 = exp{i*(M£)+^ (g(r)dY(r)-±g(rfdr 

Taking into account that 

exp {t* (r , /)} = e^ 2 / 2 e a *^e a ^\ 

we can obtain 

where f t (r,i£) = g(r) - is(r) r ,( t ), / t *(r,if) = g(r) + is(r) r ,Q) for r) t = i£. The 
integral over the trajectories of the equation (|4.9|l 

4(0 = jf-K,^ ) (t > <0)dr+ l/t|2( f ~ 9 \ 

where 

g 2 = [ 5 (r) 2 dr, |/ t | 2 ^) = [ f*(r, if) /(r, if) dr, 
J o Jo 
can be written using 1)4. 4|l in the following form 

(5.1) 7 9 (C)= / {-K+ImC^rJ.^+ierfs^s^Jdr. 

Given that e a ^'H $ = we deduce the following result. 



) 
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Proposition 3 The mathematical expectation (ip<g)6<D\Xg(t, £) tp®8<^) with respect 
to the state vector h = ip (8 5$ for arbitrary ip £ F can be written in the form 
{tp\Xg(0 V^j where the expectation xtiO °f the operator X g (t,£) with respect to 
the vacuum vector 5® € T has the form 

(5-2) X ;(C)=exp{4(e) + i?(4 9) (t,zO)}. 

The function 1 1— > Xg(£) satisfies the operator equation 

(5.3) = {(K t B4,fl)+2Im(C t fl(t),ifl+i^)}4(0, 
in the partial (functional) derivatives of the first order 

(C,d)X(0 = {i/2)(X(Z)R(Q + R(QX(€)). 

Indeed, it is not difficult to find, by means of direct substitution, that the char- 
acteristic operator- function l|5.2|l satisfies the equation (|5.3|l . However, this can be 
skipped, as far as this equation with the commutation relation 

[R(Q,X(£)]=s(C,t)X(€), vc^es, 

is concerned, it can be written in the form of the equation (|3.4[) . that was derived 
before for the vacuum expectation n g (X) of any operator of the form X g (t) = 
e g (t)X(t): 

^ M* (0 + (7 t (*(£))) = 9(t) »l(X(0 R(Q + R(Ct) X(0), 

where 

lt {X) = i [X,v t (R) + 1 -u t {R,R)] + \ (R(CI)[R(Q,X] + [X,R(d)]R(Q). 

In particular, for g = this equation defines the Markovian map /1q : A — > A for 
the quantum diffusion R(t, rj); this map also has the characteristic operator- valued 
function Xo(£) — MoC^XO) m the form l|5.2|) which is defined by the integral /q(£) 
and ^ o) (t,*f)=*0 r (t)f. 

An a posteriori linear quantum diffusion under the continuous measurement of 
the process (14 .6f) is described by the a posteriori characteristic operator-valued 
function x*(£) = M (X(£)). Its Wick symbol is the operator- valued function (|5.2(l 
defining the characteristic function 

^(0 = (V"8> /® I X (0 (V> ® /»)) exp { - Il/H 2 }, V G H, 

of g(i) = 2Kef(t) on the exponential normalised vectors /® exp{— 1 1|/|| 2 }- Con- 
sequently, the stochastic operator- function x (£) > is obtained (using a factorial sub- 
stitution in (|5.2p . with the standard Wiener process Wt replacing J Q g(r)dr ) along 
with the solution rj r = <p r {t, z£) of the backward linear stochastic equation 

(5.4) —d-Vr + iK r srj r dr = (( r + Cf) dw r , rj t = i£. 

Substituting r\ r — (f>l 9 \t,it;), defines the solution of the filtering equation for the 
case of linear diffusion and X = X(£). 
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In order to find the operator x (0 in the form of a stochastic function w ) = 
w(x (£)) (of the trajectories wt of the observable process Ij4.6|l ). let us solve the 
equation Q3.6|l for the initial Gaussian state l|4.2|) with an arbitrary #(77) = (77, #) 
instead of do(n). In that way we shall find the Wick symbol 

?'(<?,£) = K I X (0^), * - ^0 +2Re(gr7), 
where ip„ = rj® exp{— 1 ||?7|| 2 }. This symbol defines with the normal ordered 

substitution of R(r]) — do{i]) I + 2dtA(n) into 9 ($,£) instead of $(77). It can be 
obtained by solving the linear stochastic differential equation in partial (functional) 
derivatives of the first order 

idtf(0 + {s(v t7 + ^e t ( S Si,st)}e\0dt 

(5.5) = ((K tS t,d)dt + 2lm(( t dm,id+^s0)e\0, 

representing the filtering equation in terms of the nonnormalised a posteriori char- 
acteristic functional 9 (£) = <^(x*(C)) ; for the initial Gaussian characteristic func- 
tional 



Theorem 3. The solution to the filtering equation \5.5\) with the Gaussian 9 (£) 
#(£) defines the stochastici characteristic functional 9 (ui) = ui(9 ) in the form 

(5.6) o t (0 = p t ™p{$M)-\pt(tO 

In this equation = 9 (0) is the Gaussian probability density of the output process 
l4-(>\ ) with respect to Wt ■ 

f = exp I J (d r (2 Re C r ) dw r - i d r (2 Re C r ) 2 dr^ | , 

$t{v) = (Vt^t) is a linear stochastic functional of the a posteriori mathematical 
expectation of the operators R(t, n) satisfying the quantum linear filtering equation 

(5.7) dd t (n) + MiKtsn) dt = 2 Re(77, k*Cf ) dw t + v t (sr)) dt 

with the initial i?o = $ € ReS, dwt = dwt — (2Re£ t ) dt, k t = p t + (i/2)s . Here 
Pt(^v) = (tiiPtv) * s a symmetric quadratic form of the a posteriori covariance 
R(t 7 n) and R(t,£), satisfying the Riccati equation with the initial po^^rj) = (£,77) : 

d I a I 2 

(5.8) —p t (ri,r])-2pt(i],iK t T])=et{sri,sf])-\2Re(i],k t C)\ . 

Proof Let us find from l|5.5|l the stochastic equation for d (i£, uj) = In 9 (£, us) 
using the logarithmic Ito formula 

d^(i0 = J-d9\0- Udd'mf, 
9 (0 2 
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where 

(d^w) 2 = (d^(%)AX^)) 2 

= (<2ReC t ,da*(77)) + (ImC B ,S7?)) 2 dt 
= { X t (do* (r?), do 1 (»?)) - 2z P t (do* (r?) , st?) + u(sri, sr?) } dt, 
because dw t dr/3 t = dt. Here 

A t (r?,r?) = (2ReC t ,??) 2 , ^t(sry, sr?) = (ImC t , sr?) 2 , 
ftfaiO = (2ReC t ,r?)(ImC t ,0 = (iP&rj) 

This gives a quasilinear stochastic equation of the first order for a (77): 
dcr* + {s(r?,u t ) + (iK t sr?, 9?*) - -£t(sr?,sr?)} d* 

= { (2 Re C tI da ( ) + (Im (f , sr?) } dw t - ^ X t (dd\ 5a 4 ) dt, 

where Kj = Kt — (3 t and £4 = e± — v t . This equation has a solution in quadratic 
form 

o-*(r?) = In p t + d t (ri) + - p t (r?,r?), 9cr*(r?) = tf t + p t r?, dV = p t , 
where 

dmp* = da t (0) = (2ReC t ,^)dw t -iA t (S t ,^)^, 

d^ = d9a t (0) = (2p t RcC f +sImC t )dw t -{(p t A t -isKj)^_ SU( }dt 
dp t = d9 2 CT* = -{p t A t p t + se t s + i (p t K t s - sK t T p t )} dt, 

and K T (pr?,sr?) = rt(sr?, pr?). By explicit expression of the finite-dimensional op- 
erators j3 t , At and v% in terms of Ct>Ct > the stochastic integral lnp* = L d<7*(0) 
can be represented in the form given by theorem 3. The differentials d<9<7'(0) and 
dd 2 o t can be written in the form of the stochastic differential equation (|5.7|l for 
i?t(r?) = (77, -d t ) , and in the form of the Riccati equation (|5.8|) for pt(r], r?) = (r?, p*r?) 
, with k t = p t + (i/2) s because s T = — s . Hence 

(r?,2 Pt ReC t + slmQ = 2Re(r?, (p* + l - s) C\), r?eRe~°, 
h(PtV> PtV) ~ 2i/3 t (p t r?, sr?) + i/ t (sr?, sr?) = (2 Re(r?, k t C?)) 2 . 
Theorem 3 is proved. 

Remark 3 The quantum linear filtering equation 1)5. 7|) , (|5.8|l can be written in the 
form of the classical Kalman-Bucy filter 

d^M + Kr9tM+su f )dt = 2Re(kfCf)dw t! 

(5.9) — p t + a t T p t + p t a t + se t s = p f AtPt, 
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where at = XtPt + i Kts, a] = ptXt — isn 1 , Xt = 4Re£ ( Re£j . The solution of 



this system of equations with the initial $o (u) — po = 1 gives a (rj) in the form 
of the integral 

(5.10) + i [£ r (sry r ,s77 r ) - A r (p r (?7 r ),Pr (%.))] dr}, 

over the stochastic trajectories 7y r = (j> r (t,rj) of the adjoint equation l|5.4fl with 
p r = $ r + p r rj r , rj = <j> (t, rj), and 

dT t = v t dt + lm(ldw t - 

Indeed, if d-rj r = rj r — rj, dr is the backward stochastic differential, then d(i? r , rj r ) — 

(rj r , d# r ) + (i? r ,d_ry r ), and 

d [Pr(f? r ,?7 r )] = 2p r (r? r! d_r? r ) + p r {T] r ,rj r ) dr. 

Using the equations l|5.9|) and writing the equation l|5.4|) in the form d-?j r = 
a r rj r dr — 2 Re £ r du r with respect to 

du r = dw r + (p r rj r , 2 Re C r ) dr, 

one can obtain, by integrating by parts, the difference ct* (77) — hip* — i3(t] ) — \rj^: 

1 1 2 



I (77, d#) + - pr/dr - (dJfj, p(rj)) } 
= J I (s?7, dT) - (2 Re Cdw, 3) + (a??, pr?) + i (?7, rj) - A(pr), pr?) drj 
= jf { (sfj, dT) - 0(2 Re C) dw + - [e(srj, sfj) + X0, #) - X(p, p)] } 

= J {(s^,dT) + i [e(sJ7,s»7) - X(p,p)] drj - lnp*, 

which gives {5~TU|) with lnp* = J* 0(2 Re C) dui - \ X0, d) dr). 
Example. Let us consider the quantum linear filtering equations (|5.9|) for the case 
of indirect nondemolition observation of the coordinate of free quantum Brownian 
motion, described in the above example. Since £ t = Q we obtain (3 t — 0, j t = 0, 

v>t = 0, and At = A ^ 1 / = This allows us to define the unique positive 

V Reut V Imw * V 2 / 2 V ty ' 
solution of the Riccati equation l|5.9|l with the initial p = 1. Here 

, _ 1 (2 M 2 iu* t \ . f _ 1 (2 \to t \ 2 -iu> t ' 



k* - 

Recot V ~ iuj t V 2 7 ' * Rew t V iw * V 2 
correspond to the degenerate form kt (77" , rj) = \2(Jtr) p +ir] q \ 2 / (2 Reui t ) of the solution 

k t + k t Atk t = ils Q A* + k t + k t ( i A, /w, j s 
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with the initial ko = 1 + (i/2)s and fixed k( — k' = is. In addition, the parameter 
u>t € C, that defines the a posteriori wave function 



ip (x) = (v/Rewt/Trp*) cxp <{ - - u t (x - q t ) 2 + -p t x 



. t ~t 



for <p = (tp \X(£) ip ) satisfies the one-dimensional complex Riccati equation 
d i 2 1 , 1 

OJf H Ult = — A, Wn = - , 

di 2m ' 2 ' 2 
which has (only) one positive solution 

a l + ath(Af/a) /r — 777 , 
w * = 77 TTTTTTTT — ' a = ^Xm 2(1- 1). 
2 th(Xt/a) + a 

The a posteriori mathematical expectations 

1 f ~t 2 d -~t 1 — t — i 

P* = — / ^ W ~ ( x ) dx = — (V> I -FV )> 

p j z ax p 

$t = ^ / 4> t (x)*xtp t (x)dx = I QV> ), 

p .7 P 

of the momentum and the coordinate of the quantum Brownian particle satisfy the 
linear system of stochastic equations 

dp t + v A dw t = 0, p = p, 

Recut 

dq t - —p t dt = ^VA — ^ — dw t , $0 = q- 

m 2 Re cut 

The solution of this system gives the probability density 



9 



exp | J (V\q r dw r — — g?dr 



for the indirect nondemolition observation Y(t) of the coordinate of the free quan- 
tum Brownian particle with the mass m. Note, that the solution u) t tends expo- 
nentially to the limit Woo — a which defines the finite a posteriori dispersions 

^ = \\(p-pI)$ \\ 2 _ 2 K| 2 
Ho\\ 2 Rewt 

2 ||(Q-?/)^ || 2 1 



*" 2 2Rew t 



Ho 

and the correlation p t — \\ip \\ — 2 Re((P — pi) ip \(Q — ql) tp ) in the limit t — > 00: 



°lo = V2Xm, = y/2/(\m), p x = -1. 

It is well known that the unobserved free quantum particle becomes 'fuzzy' and 
has an infinite a posteriori limit dispersions. The contradiction of the localised 
motion of the quantum particle, which is seen during the continuous observation, 
and of the impossibility of localising this particle on the basis of the von Neumann or 
Lindblad equation, has provided a source for quantum paradoxes such as the Zeno 
paradox 0]. The derived quantum filtering H8.f)JI and a posteriori von Neumann 
(|3.10|l equations resolves these quantum paradoxes not only on the qualitative level, 
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but also on the quantitative level of the microscopic quantum stochastic model for 
the continuous observation. 



6. Appendix 

1. Let {Hj|i^ > 1} be a continuous family of Hilbcrt subspaces H^ C H with 
nondecreasing norms: r\ < £ HV'llr; < IIV'II?) v ' ^ G H. It will be called a scale of 
the Hilbert space H with a scalar product (ip\ip) = lim^i H^Hf. An inductive limit 
lim^ii of the Hilbert scale {H^} is defined as a pre-Hilbert space D = U^>iHj, 
provided with the inductive convergence: 

V„^0^3£> 1 : |K|| c -»0. 

The operator X : D — > D is called (inductively) continuous if — > for any 
convergent sequence {V'n} that tends to zero. This means, that the restriction of 
X to any subspace H^ € {H^} is a continuous map into some Hj C D, i.e. for any 
C > 1 there exists £ > 1, such that 

(AT) ||X||f:= S up{irak}<oo. 

The set B(D) of all continuous operators X , having the Hermitian-conjugate oper- 
ators X* on the pre-Hilbert space D form an associative algebra with the identity 
/ € B(D) and the involution X** = X. This algebra is a C*-algebra of the bounded 
operators only if D = H. 

2. Let L + be the set of all the chains T n = (ti, . . . ,t n ), ti G R+, t\ < ■ ■ ■ < t n 
of length n < oo, identified with finite subsets r C R+, t — {t\, . . . ,t n } of the 
cardinality |t| — n E {0, 1, . . .}. We denote by dr = Ylter ^ ^ ne positive a-finitc 
measure on L + = ^^Lo^n; which is defined as a sum X^^Lo - 1 "" °^ the measures 
dr„ = dti ■ ■ ■ dt n on r„ = {|r| = n} with the only atom dro = 1 on the empty 
chain r = corresponding to |r| = 0. The Hilbert space H = L 2 (T + ) of the 
square-integrable functions / : L + — > C, (/|/) < oo, where (/|/) = / |/(r)| 2 dr, 



(A.2) y s(r) dr := £ / ' ' ' / •••>*») dt i 



■<t„<oo 



is naturally identified with the Fock space T = T(JC) , of the sequences / = {VrJ n = 
0, 1, . . .} of symmetric continuations (p n : R" — > C of the functions /(r n ) , with 
the scalar product 



00 1 /■ oo /• c 

(/i/) = E^/ •••/ 

n=0 »• Jo Jo 



I 1 2 

\<p(ti, ■ ■ ■ ,t n )\ dti ■ ■ -dt n . 



The Hilbert scale > 1} of the dense subspaces 

Ft = If *T I ll/llf = / e |T| |/(r)| 2 dr < oo} 

is called Fock 0-53 over the Hilbert space K, = L 2 (R + ). The function /(r) 
5${t) ( where 5q(t) = 1 for r = 0, 8q(t) = for r = 0), normalised (||/|| 2 
/ C' T '^0( T ) dr = 1) with respect to any £, is called the vacuum function. 
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3. The most important examples of the continuous operators, and their conjugates 
T+ — > T+ on the inductive limit T+ = U^>iJF^, are the 'quantum annihilation' 
operators a t on the interval [0, t), which act as the integrals 

[2t/]W= / /(rUr)dr, V/e^. 
Jo 

Here r U r is a chain, defined almost everywhere (r ^ r) as a union {r, r} of the 
single-point chain r € R + with some chain r = {ii, . . . , t n }. Although the operators 
a t possess a complete system of exponential eigenfunctions 

f(r) = k®(T):=l[k(t), kel 2 (R+), 

t€T 

they are not normal, i.e. they do not commute with their adjoint operators ajf, 
defined by the finite sums 

r<t 

[%f\(T) = £f{T\r), V/e^+, 

where r\r = {t e r|i 7^ r}. Polynomials of the operators {a^o^i € R+} form the 
Weyl algebra over the simple functions g G L 2 (R + ), generating the whole algebra 
S(JF + ) through the triviality of the commutant: 

{Y e I X g {a*, a t } <s> = 0} = CI. 

Linear elements of the algebra B(.F+), defined as the quantum Wiener integrals 

/■OO 

- / {f(t)da t + g(t)da* t ) =a(f) +a*(g) 
Jo 

of the square-integrable functions f,g : R + — > C, together with the identity opera- 
tor 1 , form the ^representation f(/, 3)* = r(g* , f*) of the canonical commutation 
relations 

(A.3) [f(/,.g),f(/, 5 r]=(||/|| 2 -||.g|| 2 )i 

with respect to the involution (/, 3)" = (<?*,/*)• The vacuum function 5$ E J 7 
induces a Gaussian state on B{T + ), defined by the characteristic functional 

(^|e^)5 )=e xp {i(||/|| 2 -|| 5 || 2 )} 

on the Hilbert space of pairs r\ = (f,g) with the norm 2 \r]\ 2 = i (||/|| 2 + ||ff|| 2 ) 
and symplectic bilinear form s(t]' \rj) which equals s{rfi,rj) = i (||/|| 2 — ||g|| 2 ) for 

v' = (9*,f*)- 

4. The measurable function F : t 1— > F(t) whose values are the operators T> — ► T> of 
the inductive limit I? = Uj>i7i^ is called locally p-integrable in the inductive scale 
{H^ = Hj ® J^} , if, for any £ > 1, there exists a £ > 1, such that 

:= (imr)||f) P dr) ^ < ex., Vie R + . 

In particular, this condition means that the operators F(t) :T> — > T> are continuous 
for almost all t € R + . For such square-integrable functions the quantum stochastic 
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integrals a% (F), a t (F) on {a*\r < t} and {a r \r < t} are denned to be the operators 

r<t 

[a* t (F)h](r) = £[F(r)/i](-r), Vh E T> 1 
(A.4) [a t (F)h](r) = f [F(r) h(r)] (r). 

In this equation h <— > h(r) is a Maliven derivative, defined in the Fock representation 
h : T + — ► D of the elements h E T> almost everywhere (r ^ r) by the vector-function 
h(r) E T> as h(r,r) = h(r U r). The continuity of the operators at(F), a,f(F) in I? 
follows directly from the estimates 

11^(^)11^ < n/V^II!;?, Ve>0, 

||S*(F)||^ £ < VU~e\\F\\f 2 t , V £ <e 

obtained in 0, From these estimates, if there exists £ > 1 for any £ > 1 

such that H-Fll^'j, ||-D|||'? < oo then, for any ( + > 1 , there exists a £_ > 1, for 
which the operator 

i t (F,D)=a t (F)+a* t (D) 
is bounded from to . To be precise, by choosing for every £ + > 1 a £ > 1, 
such that C < C + , and an £ > £ — 1, such that e < £ + — £, we obtain 

||4(F,D)||^< + pill;?) 

for any £_ on the nonempty interval (1, £ — e], 

5. If D(t) £ B(V) for almost all t € R + and the adjoint function £>*(*) = D(t)* is 
also locally integrable on D, then a1{D) E B{V) and a£ (D)* = a t (D*). This means 
that the integral z (D*,D) is (formally) self-adjoint. Moreover, the set of integrals 

X(t) = i t (D*,D)+ [ G(r)dr, 
Jo 

where G : R+ — > £>(2?) is locally integrable (p = 1), together with the G* , func- 
tion, forms a ^-algebra with respect to the pointwise operator product (X*X)(t) = 
X(t)*X(t). This product is defined by the quantum nonadapted Ito formula 10 

X(t)*X(t) = i l {D*X + X*D,F*D + X*D) 

(A.5) + [ [G*X + X*D + D*D + D*X r + X*G](r) dr, 

Jo 

which corresponds to the case of the locally square-integrable operator function 
t h- > X t (t) 1 X t *(t) = X(t)l. Here IhI, means a derivation, which is defined by 

[X t h][T) = [Xh] (rLAt)-[Xh(t)](r) 

for almost all r E T, t ^ r. In this instance, the case X t (t) = = X£(t) for all 
i S R+ corresponds to the adaptive property of the operator functions F(t),D(t) 
and G(t) . 

Comprehensive information on explicit stochastic integration in the Fock scale 
is given in |2()j . 
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